General non-Markovian dynamics of open quantum systems 
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We present a general theory of non-Markovian dynamics for open quantum systems. We explore 
the non-Markovian dynamics by connecting the exact master equations with the non-equilibirum 
Green functions. Environmental back-actions are fully taken into account. The non-Markovian 
dynamics consists of non-exponential decays and dissipationless oscillations. Non-exponential decays 
are induced by the discontinuity in the imaginary part of the self-energy corrections. Dissipationless 
oscillations arise from band gaps or the finite band structure of spectral densities. The exact analytic 
solutions for various non-Markovian environments show that the non-Markovian dynamics can be 
largely understood from the environmental-modified spectra of the open systems. 
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Understanding the dynamics of open quantum systems 
is one of the most challenging topics in physics, chemistry, 
and biology. The environment-induced quantum dissipa- 
tion and dccoherence dynamics are the main concerns in 
the study of open quantum systems 0, 0j ■ Decoherence 
control has also recently become a key task for practical 
implementations of nanoscale solid-state quantum infor- 
mation processing Q, where the decoherence is mainly 
dominated by non-Markovian dynamics due to the strong 
back- actions from the environment. A fundamental issue 
is how to accurately take into account non-Markovian 
memory effects, which have attracted considerable at- 
tention very recently both in theory (34l2j and in exper- 
iments [13l4l5| . 

The non-Markovian dynamics of an open quantum sys- 
tem can be described by the master equation of the re- 
duced density matrix p{t). This is obtained by trac- 
ing over the environmental degrees of freedom, p(t) = 
tr[ptot(i)]j where ptot(^) is the density matrix of the total 
system. The standard approach to the non-Markovian 
dynamics uses the Nakajima-Zwanzig operator projec- 
tive technique [l6| where the master equation is formally 
written as 



dp® 
dt 



drK.it - r)p(r). 



(1) 



The non-Markovian memory effects are taken into ac- 
count by the time non-local integral kernel K.{t — r). In 
practice, very few systems can be exactly solved from ([1]). 
Therefore, the generality of non-Markovian dynamics has 
not been fully understood. 

In general, there are three typical timescales in an open 
system to characterize non-Markovian dynamics: (i) the 
timcscale of the system ~ l/e s , where s s is a typical 
energy scale of that system; (ii) the timcscale of the en- 
vironment 1/d, where d is the bandwidth of the en- 
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vironmental spectral density; (iii) the mutual timescale 
arising from the coupling between the system and the 
environment ~ where T is the dominant coupling 
strength. It is usually believed that non-Markovian mem- 
ory effects strongly rely on the relations among these dif- 
ferent timescales. However, such relationships have not 
been quantitatively established yet. 

Here, we show that the general behavior of non- 
Markovain dynamics is mainly determined by the band 
structure of the environment and the coupling between 
the system and the environment. We explore the non- 
Markovian dynamics from the analytical solution, solved 
by connecting the exact master equation with the non- 
equilibirum Green functions. Exact master equations 
have been derived only for limited system-environment 
couplings @, QjHlII]- Establishing the connection be- 
tween the master equation and the non-equilibrium 
Green functions provides a new way to explore the non- 
Markovian dynamics even if the exact master equation 
of the open system is unknown. 

Exact master equation and non- equilibrium Green 
functions. — We begin with a fermonic (bosonic) many- 
body system consisting of N single-particle energy lev- 
els £, (i = 1,2,--- , N), coupled, via particle-particle ex- 
changes, to a fermonic (bosonic) environment, H$b = 
J2akii V akialb a k + V* ki b\ k ai}. Tnc environment can con- 
tain many different reservoirs, each reservoir is specified 
by its spectral density J alJ (uj) = ^J2k V akiV* kj S(uj- 
efe), where V a H is a coupling strength between the system 
and reservoir a. The operators a\ (of) and b* ak (b a k) are 
the particle creation (annihilation) operators of the di- 
crete energy level i of the system and the continuous level 
k of reservoir a, respectively. These creation-annihilation 
operators obey the standard anticommutation (commu- 
tation) relationship for fermions (bosons). Nonlinear 
particle-particle interactions in the system are not in- 
cluded. Using the coherent-state path-integral method 
[20| to the Feynman- Vernon influence functional [2l[ , the 
exact master equation of such an open system can be de- 
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dp® 
dt 



1 r 

i 



Hs(t),p(t) 



2a,jp(t)a 



ij 

- a\ajp{t) - p{t)a\a^ + j t] (t) a\p{t)aj 
±a j p{t)a\-a\a 3 p(t)Tp{t)a ] a\ }. (2) 



The first term in © is the unitary term with the renor- 
malized Hamiltonian Hs(t) = ^2ij£sij(t)a\aj. The sec- 
ond and third terms give the non-unitary dissipation and 
fluctuations, respectively. The ± and =F signs in the third 
term correspond to the system being bosonic/fermionic. 
The renormalizcd energy levels s s (t) , the time-dependent 
dissipation coefficient ~f(t) and the fluctuation coefficient 
7(i) in © are given by 



e s (t) = -[u(t,i )-it- 1 (i,i )-H.c], 



7 (i) =--[u(t,*o)u- 1 (t,*o)+H.c.] ) 



(3a) 
(3b) 



7^*) = »(V)-[*(*,to)u -1 (*,toMM) + H.c.]. (3c) 

In Eqs. ([3]), the N x N matrix functions u(t, to) 
and v(t,t) are related to the non-equilibrium Green 
functions of the system in the Schwinger-Keldysh non- 
equilibrium theory [H, H3|, u io (t, t ) = ([^(i), a](f )] T ), 
&ndvij(t, t) = {a]j(i)ai{t)) subtracting an initial-state de- 
pendent part [iH . These Green functions obey the Dyson 
equations, 



dr 



u(t, t ) + i£ s u(t, t ) + / dr'g(T, t')u(t\ t ) = 0, 



to 



(4a) 



dr 



v(t, t) + ie s v{r, t) + / dr'g{T, t')v(t', t) 



to 



[dr'g(T,T')u^T',t ), (4b) 

Jtn 



'to 

subjected to the boundary conditions u(to,to) = 1 and 
v(to,t) = with to < t < t, where £ s is a N x N ma- 
trix given by the bare single-particle energy levels of the 
system. The self-energy corrections, g(r, r') and g(r, r'), 
which take into account all the back-actions from the en- 
vironment, are expressed explicitly by 



/ |£j a (a,)e-M— ') , ( 5 a) 



(5b) 



known exact master equation and the non-equilibrium 
Green functions for open quantum systems. 

General non-Markovian dynamics. — Different from 
the Nakajima-Zwanzig master equation, the exact master 
equation ([2]) is local in time, characterized by the dissipa- 
tion and the fluctuation coefficients, j(t) and 7(f). Non- 
Markovian memory effects are manifested as follows: 

(i) The coefficients j(t) and 7(t) arc microscopi- 
cally and non-perturbatively determined by the non- 
equilibrium Green functions from the Dyson equations 
(|4|) . The non-Markovian memory effect is fully coded into 
the homogenous non-local time integrals in Q with the 
integral kernel g(r, t'). In other words, the self-energy 
correction g(r, t') serves as a memory kernel that count 
all the back-actions from the environment. 

(ii) The coefficients 7(4) and ■y(t) are constrained 
by the non-equilibrium fluctuation-dissipation theorem. 
The inhomogenous non-local time integral in (|4b[) with 
the integral kernel g(r, r'), depicts the fluctuation arisen 
from the environment. Because v(to, t) = 0, we can ana- 
lytically solve Eq. (|4b|> : 



v(r,t) = / dr\ I dT 2 u(T, n) g(n, t 2 ) M f (t, r 2 ) 

J t J to 



(0) 



This solution shows that Eq. (|3cj) is a generalized non- 
equilibrium fluctuation-dissipation theorem in the time 
domain (the reduction to the equilibrium fluctuation- 
dissipation theorem is given in [251]). The fluctuation- 
dissipation theorem is a consequence of the unitarity of 
the whole system. It guarantees the positivity of the 
reduced density matrix during the non-Markovian time 
evolution. 

Based on the above intrinsic features of open quantum 
systems, we can now explore the general properties of 
non-Markovian dynamics. From Eqs. ([3]), we can express 
the Green function u(t, to) in terms of the dissipation 
coefficient -f(t) as 



u{t,t ) 



Texp |-jTdr[ie(r) + 7(7-)] }, 



(7) 



where T is the time-ordering operator. This solution 
indicates that u(t, to) fully determines the dissipation 
dynamics of the system. However, due to the time- 
dependence of the dissipation coefficients, the detailed 
dissipation dynamics can vary significantly for different 
environments. 

Explicitly, equation (JSJ show that g(r, t') = g(r— r') 
and g(r, t') = g(r — r'). Thus we can write u(t, to) — 
u (t — to). Using the modified Laplace transform U(z) — 
J^dtu(t)e lz ^~ to \ it is easy to obtain 



where the function f a (ui) 



T I]" 



is 



the Bose-Einstein (Fermi-Dirac) distribution of bosonic 
(fermionic) reservoir a at the initial time t Q . Equa- 
tions ©-© establish a rigorous connection between the 



U(z) = — 
zl ■ 



(8) 



e. -£(*)' 

where I is the identity, is the Laplace transform of 
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FIG. 1: (color online) A schematic pole structure of the Green 
function U(z). The thick red lines on the real z axis corre- 
spond tO J2a Ja{z) =fi 0. 



the self-energy correction, 

did J a (uj) z=U j±io+ 



2tt z — uj 



A( W )Tf£ 



J a (uj) 



(9) 



and A(w) = Ea^J If 4^ is the Principal value of 
the integral. It can be shown that the general solution of 
u(t,to) is given by 

u(t - to) = E ^ie-^-^ + E / ^ f* 7 ^ + i0+ ) 



C/(o;-i0 + )le- iw (*- to ). 



(10) 



The first term in (|10|) corresponds to localized modes 
with poles {w,;} located at the real z axis with 
J2 a Ja{u) = 0. The coefficients {Zi} are the corre- 
sponding residues. The localized modes exist only when 
the environmental spectral density has band gaps or a 
finite band, i.e., ^2, a Ja(w) vanishes in some frequency 
regions, see Fig. [TJ These localized modes do not decay, 
and give dissipationless non-Markovian dynamics. The 
second term in (|10[) is the contribution from the branch 
cuts {-Bfe}, due to the discontinuity of E(z), so does U(z), 
across the real axis on the complex space z, see Eq. (|9|). 
The branch cuts usually generate non-exponential decays 
[27| . which is another significance of the non-Markovian 
dynamics. When the system is weakly coupled to the 
environment, the non-exponential decays are reduced to 
exponential-like decays. 

Equation (fTU)) provides indeed a general solution of 
the non-Markovian dissipation dynamics. It shows that 
the non-Markovian dissipation dynamics consists of non- 
exponential decays plus dissipationless localized modes. 
Such a solution for the two-point Green function u(t,to) 
is generic and can be proven from the quantum field 
theory [33| . even if particle-particle interactions are in- 
cluded. 

The Green function u(t,to) reveals the general non- 
Markovian dissipation dynamics. The non-Markovian 
fluctuation dynamics is constrained by the fluctuation- 
dissipation theorem via the Green function v(t,t) of 



(JH). Thus, the whole picture of non-Markovian dynam- 
ics is fully characterized by the dissipation and fluctua- 
tion coefficients of ([3]). The non-exponential decay part 
of pop makes the dissipation coefficient ~f(t) oscillates 
between positive and negative values, representing the 
back-flow of information from the system to the envi- 
ronment [1, H|. Non-exponential decays alone give -y(t) a 
non-zero asymptotical value. If there are localized modes, 
•y(f) will vanish in the steady-state limit, resulting in 
dissipationless oscillations. In the weak coupling region, 
"fit) can still be time-dependent but keeps positive, the 
corresponding dynamics gives simple exponential-like de- 
cays, observed mainly in the Markovian limit. Further- 
more Eqs. (|3"c|) and (JB]) together show that except for 
the initial environmental temperature dependence, the 
time-dependence of the fluctuation coefficient j(t) be- 
haves similar to fit), due to the fluctuation-dissipation 
theorem. In conclusion, non-Markovian dynamics can be 
fully understood from the solution of the Green function 
u(t,t ). 

Examples and discussion — To be more specific, let us 
first examine the non-Markovian dynamics of a single- 
mode bosonic nanosystem, such as a nanophotonic or 
optomechanical resonator, coupled to a general non- 
Markovian environment with spectral density 



J(cj) = 2irr]uj[ — 



u> 

exp ( 

UJ C 



(11) 



where r\ is the coupling constant between the system and 
the environment, and uj c is the frequency cutoff. When 
s = 1, < 1 and > 1, the corresponding environments are 
Ohmic, sub-Ohmic and super-Ohmic, respectively [26| . 
Following the above general procedure, the analytical 
solution of the non-Markovian dissipation dynamics is 
given by (setting ip = for simplicity): 



u(t) = Ze'*"'* + - 

7T 



duj- 



J(i>j)e~ 



A\uj - e s - A(uj)} 2 + ,P{uj) ' 

(12) 



where A(w) = ±[£(w+i0 + ) + X;(w-i0 + )] and the Laplace 
transform of the self-energy correction 



r]ijj c \jv^/—uje "erfc(\/— uj) — t/it\ s = 1/2 
r\u) c [to exp(— w)Ei(w) — l] ,s = 1 

^ r]uj c [£ 3 e- s Ei(w) - uj 2 - uj - 2] s = 3 



(13) 



with uj = lo/uj c . Due to the vanishing spectral density for 
u) < 0, a localized mode at uj' — e s — £(u;') < occurs 
when ijuj c T(s) > e s , here T(s) is a gamma function. The 
localized mode leads to the dissipationless process. The 
corresponding residue is Z = [1 — E'(a;')] _1 . This ana- 
lytical solution precisely reproduces the exact numerical 
solution in the previous work [Io[ • Figure [5] shows that 
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for a small 77, the dissipation dynamics is an exponential- 
like decay, The corresponding j(t) and j(t) are time- 
dependent but positive (corresponding to Markovian dy- 
namics). When 77 > 0.3, the non-exponential decay dom- 
inates, and 7(7;) and 7(f) oscillate in positive and negative 
values with nonzero asymptotical values. When 77 > 0.6, 
the localized state occurs, and u(t) does not decay to 
zero. Correspondingly, 7(f) and 7(7;) asymptotically ap- 
proach to zero. 
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FIG. 2: (color online) The time evolution of the Green func- 
tion u(t), the dissipation and the fluctuation coefficients, 7(7;) 
and 7(i), in a sub-Ohmic bath, for several different val- 
ues of the coupling constant 77. We take other parameters 
e s = 13.83/jeV, uj c — e a and also fcsT = e a . 

The second example is a fermonic system, a single 
electron transistor in nanostructures which consists of 
a quantum dot coupled to a source and a drain. The 
source and the drain are treated as two reservoirs of the 
environment. Their spectral densities take a Lorentzian 
form with a sharp cutoff, 



r d 2 

a 11 a 



(w - ojc) 2 + d 2 a 



Q(Q — |cj — uj c 



(14) 



with a = L(R) for the source (drain), where d a is the 
halfwidth of the spectral density and T a is the coupling 
strength between the system and reservoir a. We add 
a sharp cutoff to simulate a finite band for the environ- 
mental density of states. When CI — > 00, the above spec- 
tral density is reduced to the usual Lorentzian spectral 
density that has been used in various studies of nano- 
electronics fA [28l - [3l| . We consider the symmetric case, 
(rj,,dz,) = (TR,dii) = (T,d). Then the analytical solu- 
tion of the Green function u(t) becomes 



*(*) 



2 
i=i 



IOJa t 



1 



u> +Q 

du- 



J(ui)e 



A(w)] 2 + J 2 (w)' 
(15) 



where J(u>) = Jl(w) = Jr{oj) and A(w) is the real part 



of the self-energy 



E(w) 



LO r — CI — LU 2(ui — UJr) 1 CI 

log -I— + ~\ tan" 1 - 

lo c + S I — u! d d 



(16) 



The two localized states are located outside of the band 
i.e., u>j = e s + T,(u)'j), with uj[ < ui c — Cl 

[1 



and u>2 > 



Cl. The corresponding residue is given by Zj = 



s'K)]- 



Again, the localized modes lead to a 



dissipationless process and the integral term shows a 
non-exponential decay. Taking Cl — > 00, the two local- 
ized modes are excluded, and the solution of u(t) repro- 
duces the exact non-Markovian dynamics of the usual 
Lorentzian spectral density (for detailed derivation, see 

The third example is a two-level system involving 
single-photon processes (spontaneous emission). In gen- 
eral, a multi-level atomic open system does not obey the 
master equation However, the Schrodinger equa- 

tion of a two-level atomic system with only spontaneous 
single-photon emission processes (at zero temperature) 
can be reduced to the Dyson equation of (|4"a|) (341436T ] . 
For a two- level artificial atom, such as a quantum dot, 
embedded in photonic crystals, because of the photonic 
band gap it was shown [351 ] that the corresponding so- 
lution contains exponential decays, non-exponential de- 
cays, and localized bound modes all together. We find 
analytically |37| that the complex pole with exponential 
decay shown in (35| has been included in the brach-cut 
integral of (fTDj) . Explicitly, the spectral density of the 
photonic crystals takes J(oj) = ^== 6(0; — ui e ) [34|,[35|. 
From Eq. (|T0|) . wc directly obtain the analytical solution 
of the spontaneous emission dynamics 



t (t) 



2uj r 



iiOr + A 



e *™). + £/2 

7T 



{Lu-e s ) 2 {Lu-Lu e ) + C 21 
(17) 



where w r is the real root given by (u> r + A)^/w7 = C, and 
A = e s — uj e is the detuning. This analytical solution re- 
covers both the exact analytical and numerical solutions 
given in [HI, Hf| . 

The above examples show that very different open sys- 
tems coupled to very different environments obey the 
same solution, Eq. (fTU|) . of the non-Markovian dynamics. 
The solutions of these examples can further be written 
in general as u(t — to) = ^T>(uj) cxp{— iuj(t — £ )} 
with 



V{uj) = 2tt^2 Zj5(u + 



J(u) 



-A(w)] 2 + J 2 H/4' 
(18) 



Equation (|18j) shows that the environment modifies the 
system spectrum as a combination of localized modes 
(dissipationless process) plus a continuum spectrum part 
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(non-exponential decays). Remarkably, the result ob- 
tained from these simple examples gives indeed the gen- 
eral structure of two-point correlation functions in arbi- 
trary complicated systems, see Ref. (33[. This indicates 
that alternatively, non-Markovian dynamics can be fully 
characterized by the environmental-modified spectrum of 
the system. If the spectrum of the open system can be 
measured, the non-Markovian dynamics can be extracted 
from its Fourier transform. This largely simplifies the ex- 
ploration of the general properties of non-Makovian dy- 
namics for more complicated open systems. 

Conclusion. — By connecting the exact master equa- 
tion with the non-equilibrium Green functions in many- 
body systems, we derive a general analytical solution 
of non-Markovian dynamics for open quantum systems, 
i.e., Eq. (fTU)) or (jT5J) . From the analytical solution, we 
show that the general non-Markovian dynamics consists 
of non-exponential decays and dissipationless oscillations. 



The dissipationless processes arise from band gaps or fi- 
nite band structures of the environmental spectral den- 
sities. The non-exponential decays are induced by the 
discotinuity in the imaginary part of the self-energy cor- 
rections from the environment. The exponential decays 
observed in Markovian limit is a special case in the weak 
coupling limit. Since the non-equilibrium Green func- 
tions are well-defined for arbitrary quantum systems, this 
theory may also provide a new approach to explore non- 
Markovian dynamics for more complicated open systems 
whose exact master equation may be unknown. 
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